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Abstract 

In this contribution we briefly review the current status of the dipole mod- 
els and parton saturation on the basis of results presented at the HERA-LHC 
workshops in the years 2006-2008. The problem of foundations of the dipole 
models is addressed within the QCD formalism. Some limitations of the mod- 
els and open problems are pointed out. Furthermore, we review and compare 
the currently used dipole models and summarise the applications to describe 
various sets of HERA data. Finally we outline some of the theoretical ap- 
proaches to the problem of multiple scattering and saturation. 

1 Introduction 

Dipole models [1-3] represent a QCD motivated framework that has been successfully applied to de- 
scribe a variety of gluon mediated scattering cross sections at high energies. In particular, they provide a 
transparent and intuitive picture of scattering processes. Their main strength is a combination of univer- 
sality, simplicity and efficiency. The dipole models are capable of simultaneously describing all F2, Fl 
and heavy quark production ep data at small x, the inclusive diffractive data, the bulk of measurements 
for exclusive diffractive vector meson production, deeply virtual Compton scattering (DVCS), and even 
nuclear shadowing [4-13]. This unified description is achieved using only a few parameters with a trans- 
parent physical meaning, such as the normalisation of the gluon distribution at a low scale, the quark 
mass or the proton size. At the same time, the dipole models provide a phenomenological insight into 
important aspects of high energy scattering, like the relative importance of multiple scattering or higher 
twist contributions. This importance may be quantified in terms of a saturation scale, Qs, the scale of 
the process at which the unitarity corrections become large [4]. Up to now, the dipole models applied to 
HERA data offer one of the most convincing arguments for the dependence of this scale on the scattering 
energy and provide one of the best quantitative estimates of the saturation scale [4-6, 11,12]. This shows 
the complementarity of dipole models to the rigorous framework of coUinear factorisation, within which 
the description of multiple scattering, although possible in principle, is quite inefficient. It is not only 
very demanding from the technical side (for instance, even the basis of twist-four operators is not fully 
understood yet), but it would also require introducing a set of new unknown functions parameterising 
the expectation values of higher twist operators at the low (input) scale. In dipole models this problem is 
bypassed by simply fitting the (implicitly) resummed multiple scattering cross section together with the 
nonperturbative contribution with constraints imposed by the unitarity of the scattering matrix. 



2 Foundations and limitations of dipole models 

Let us consider a 2 ^ 2 scattering amplitude of i+p — > f +p , where the strongly interacting projectile 
i hits a hadronic target p and undergoes a transition to a state /, while the target scatters elastically. At 
HERA the projectile is always a virtual photon, 7*, with a four-momentum q and virtuality = —Q^, 
and the target is a proton, with initial momentum p and final momentum p' . The final states considered 
are virtual and real photon states, vector meson states and diffractive states. The states i and / carry a 



Fig. 1: High energy scattering in the dipole representation. 



typical scale Q^; for i = / = 7*(Q^), = Q^. The invariant collision energy s = {p + q)"^ is assumed 
to be large, s ^ and s » where t = {p — p')^ is the momentum transfer. We shall also use the 
variable x = Q'^/s, that reduces to the Bjorken x for the case of deeply inelastic scattering (DIS). 

The key idea behind dipole models is a separation (factorisation) of a high energy scattering ampli- 
tude, A^^^'^''^, into an initial (^'j) and final (^'/) state wave function of the projectile i and the outgoing 
state /, and a (diagonal) universal scattering amplitude of a multi-parton Fock state, off a target p; see 
Fig- El The scattering operator, T, is assumed to be diagonal in the basis of states that consist of a definite 
number of partons, n, with fixed longitudinal momentum fractions, (A; = 1, . . . , n), of the projectile, 
definite helicities, Afc, and transverse positions, r^. One may write symbolically (see e.g. [14]): 

j^ip^fP = j [d^Vfe] j Kzfe] ^>){n, {zk,rk, \k}) T{J^n) ^i(n, {zfc, r,,, Afc}). (1) 

n,.7^n,{Ak} 

In most practical applications one takes into account only the lowest Fock states, composed of a quark- 
antiquark (qq) pair and, possibly, one additional gluon (qqg). In the limit of a large number of colours, 
Nc — > oo, flavourless scattering states, i and /, may be represented as a collection of colour dipoles [2]. 
For the simplest case of qq scattering, the intermediate state J^2 is defined by the quark and antiquark 
helicities, the longitudinal momentum fraction, z, of the projectile carried by the quark, the dipole vector, 
r = r2 — Ti, and the impact parameter vector, b = zri + (1 — z)r2. It is convenient to define 
the imaginary part of the dipole scattering amplitude (assuming independence of the azimuthal angles), 
J\f{x,r, b) = lm.T{T2), and the 5-dependent dipole-target cross-section 

^ = 2Mix,r,b). (2) 

The picture encoded in ([T]l may be motivated within perturbative QCD. In the high energy limit of 
QCD [15, 16], the dominant contribution to scattering amplitudes comes from vector boson (gluon) ex- 
changes, that lead to cross-sections constant with energy (modulo quantum corrections that may generate 
an additional enhancement). For each spin- 1/2 fermion (quark) exchange in the t-channel the amplitude 
is power suppressed by a factor of In consequence, the high energy scattering amplitude may be 

factorised into the amplitude describing slow (in the target frame) gluon fields and the amplitude of fast 
parton fields of the projectile moving in the gluon field of the target. This is, in fact, the basic assumption 
of the kx- (high energy) factorisation [16, 17]. In the high energy limit, the vertex describing the coupling 
of the fast s-channel parton (quark or gluon) to a gluon exchanged in the t-channel is eikonal: the large 
light-cone component of the longitudinal parton momentum and the parton helicity are conserved. Also, 
up to subleading terms in the collision energy, the fast parton does not change its transverse position in 
the scattering process. These properties of high energy amplitudes in QCD were used to derive the dipole 
model for hard processes. In more detail, the scattering amplitudes in the dipole model follow from the 
QCD scattering amplitudes obtained within the /cy-factorisation scheme, in the high energy hmit and at 
the leading logarithmic (LL) ln(l/x) approximation [1]. 
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The fact that the QCD dipole model follows from the Ajr-factorisation approximation implies that 
the model, up to subleading terms in 1/s, is also consistent with the leading order (LO) coUinear approx- 
imation [17]. In addition, as in the case of the A;r-factorisation framework, the dipole model incorporates 
an exact treatment of the quark transverse momentum in the box diagram. These kinematic effects, when 
analysed within the collinear approximations, manifest themselves as higher order corrections to the co- 
efficient functions [17]. Although the implicit resummation of the collinear higher order terms in the 
dipole model is only partial, it should still be viewed as an improvement of the LO collinear approxima- 
tion. 

Practical use of dipole models is not restricted to hard processes, where precise predictions can 
be obtained within the collinear factorisation framework. On the contrary, one of the most successful 
applications of the dipole model (the saturation model [4]) provides an efficient and simple description 
of the transition from the perturbative single scattering regime (the colour transparency regime) to the 
multiple scattering regime as a function of the process scale and scattering energy (or and x). In this 
transition region scattering amplitudes are expected to receive contributions both of the nonperturbative 
nature and from perturbatively calculable multiple scattering effects. The nonperturbative effects in high 
energy scattering are currently not computable with theoretical methods and have to be modelled. The 
multiple scattering effects enter the scattering amplitudes e.g. as higher twist contributions [isll], that 
are suppressed by inverse powers of the hard scale, Q^, and additional powers of Nevertheless, the 
higher twist effects may be quite sizable at small x and at moderately large [18]. This originates from 
a rapid growth of the multi-gluon density with decreasing x: assuming the large Nc limit, the n-gluon 
density evolves approximately as the single gluon density to power n [19, 20]. Thus, at decreasing x the 
multiple scattering effects are increasingly enhanced and may eventually become competitive with the 
single scattering contribution. 

Thus far we discussed the dipole model from the perspective of perturbative QCD. An interesting 
attempt to provide foundations of the model in a general (i.e. non-perturbative) framework was recently 
put forward [21, 22]. The scattering amplitudes were written in terms of skeleton diagrams and the 
QCD path-integral. Approximations and assumptions necessary to recover the dipole model amplitudes 
were identified. To a large extent the conclusions from that analysis confirm those obtained within the 
perturbative framework: the dipole model accuracy is not theoretically guaranteed when higher twist and 
higher order corrections are large. An interesting point raised in Refs. [22, 23] is the dependence of the 
dipole cross section, cjgg, on the dipole-target collision energy, y/s. In most models one assumes that 
Gqq depends on s through x = Q"^ /s. The scale, however, is part of the wave functions and it is not 
obvious that the dipole cross section should depend on rather than on the dipole variables, like e.g. 
the dipole scale, Interestingly, assuming the dependence of Ugg on a combined variable sr^ was 

shown to create some tension between the HERA data on F2 and Fl and the dipole model, irrespective 
of the detailed functional form of aqq. Some insight may be gained from inspecting the issue in the 
/cT-factorisation approach. Then, the energy dependence enters through Xg of the gluon, that essentially 
depends on the external state virtuality, the scattered quarks' transverse momenta and the distribution 
of the quark longitudinal momentum. So, the proposed replacement of by might be somewhat 
oversimplified. On the other hand, within the LL(l/x) approximation the standard choice of Xg ~ Q"^ /s 
is justified. To sum up, the choice of the optimal dimensionless variable that would carry the energy 
dependence of the dipole cross-section remains an open and interesting problem. 

3 Phenomenology of dipole models 

Implementations of multiple scattering in colour dipole models are based on two main approaches, 
that adopt different approximations. The Glauber-Mueller (GM) eikonal approach [24] is used in the 
family of models that evolved from the Golec-Biernat-Wiisthoff (GBW) model [4]. One assumes in 

'Multiple scattering effects that occur at low scales are absorbed into the input gluon density at the initial scale. 



3 



(a) 



(b) 




Fig. 2: The 7*p scattering amplitude with unitarisation achieved via (a) eikonal diagrams or (b) fan diagrams. For exclusive 
diffractive processes, such as vector meson production {E = V , J /tji, (j), p) or DVCS {E = 7), we have x' ^ a; <g 1 and 
t — {p " p')^- For inclusive DIS, we have _E = 7*, a; = a;' <C 1 and p = p' . 



this approach that muhiple colour dipole scatters are independent of each other, see Fig. [2^. This as- 
sumption may be supported (ahhough it was not yet exphcitly derived) with properties of the collinear 
evolution of quasi-partonic operators describing the multi-gluon density in the proton, and in the large 
Nc limit [18-20]. Assuming in addition a factorised ^-dependence of the gluon distribution, one postu- 
lates the dipole-proton scattering amplitude of the form: 

M{x, r,b) = l- exp (^-^r^asifi^) xg{x, /x^) r(6) j , (3) 

where the scale = C/r^ + //q with /io ~ 1 GeV. HERA data on exclusive vector meson production 
imply a Gaussian form of the proton shape in the transverse plane, T{h), with ylJP^ = 0.56 fm. The 
corresponding quantity determined from the proton charge radius (0.87 fm) is somewhat larger, \/lJ^ = 
0.66 fm, implying that gluons are more concentrated in the centre of the proton than quarks. The form 
^ is denoted by the "b-Sat" model [6, 11]. It can be considered to be an improvement on a previous 
model [5] where T{b) oc Q{Rp — h) was assumed, and also on the original GBW model [4] where 
additionally the scale dependence of the gluon distribution was neglected, that is, xg{x^ /z^) oc x^'^ was 
assumed for a fixed power A ^ 0.3. Note that in the GBW model large saturation effects were needed 
to get from the hard Pomeron behaviour (~ x^^'^) at small dipole sizes to soft Pomeron behaviour 
(~ x~^-^) at large dipole sizes. On the other hand, in Refs. [5,6, 11] this transition can alternatively be 
achieved with DGLAP evolution, therefore saturation effects are correspondingly smaller. 

In the alternative approach one exploits solutions of the Balitsky-Kovchegov (BK) equation [25]. 
It was derived for scattering of a small colour dipole off a large nucleus, composed of A nucleons. 
The LL BK equation rigorously resums contributions of BFKL Pomeron fan diagrams (Fig. I2J5), that 
are leading in A, 1/Nc and in the Inl/x approximation (properties of solutions of the next-to-LL BK 
equation are not known yet and so cannot be used in the dipole models). A colour glass condensate (CGC) 
dipole model parameterisation [8] was constructed from an approximate solution of the BK equation: 

„/ I ln(2/rQa) ' 



Mix,r,b)=T{b)Mix,r) = Q{Rp-b){^^[~) ' "^^^ < ^ ^ (4) 



where Qs = {xq/x)^/"^ is a saturation scaleU The original analysis [8] neglected the charm quark 
contribution to F2. The inclusion of charm was later found [11] to significantly lower the saturation 
scale when the anomalous dimension 7^ was fixed at the LO BFKL value of 0.63. By letting 7^ go 
free, a solution was subsequently found with 7^ = 0.74 which included heavy quarks but had a large 
saturation scale [9]. (This model has been modified to include a t dependence in the saturation scale 



^In what follows we shall use Qs (with a lower-case s) to denote the saturation scale defined in a model-dependent way. 
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allowing the description of exclusive diffractive processes [10].) However, the HERA data do not show 
a strong preference for the solution with 7^ = 0.74, and a secondary solution with 7^ = 0.61 and a 
much smaller saturation scale also describes the data well [12]. The CGC model dUl assumes a factorised 
b dependence which is not supported by HERA diffractive data, where one finds a significantly non- 
zero effective Pomeron slope Op, indicating correlation between the b and x dependence of the dipole 
scattering amplitude. A more realistic impact parameter dependence was included by introducing a 
Gaussian b dependence into the saturation scale Qs, denoted by the "b-CGC" model [11, 12]. It was not 
possible to obtain a good fit to HERA data with a fixed 7^ = 0.63 [11], but on freeing this parameter, 
a good fit was obtained with a value of 7^ = 0.46 [12], close to the value of 7^ ~ 0.44 obtained 
from numerical solution of the BK equation [26]. However, the value of A = 0.119 obtained from the 
"b-CGC" fit [12] is lower than the perturbatively calculated value of A ~ 0.3 [27]. 

In both the approaches to unitarisation one neglects multi-gluon correlations in the target. Thus, 
the key difference between the eikonal and the BK approaches is that in the latter one resums the leading 
logarithms of 1/x while in the former one aims at keeping a reasonable representation of leading loga- 
rithms of Q^. Both dipole model realisations have built in saturation of the black disc limit of the colour 
dipole scattering amplitude. This means that the absolute value of the T-matrix elements tends to unity 
for large dipoles or as x ^ 0. It is curious that the choice of approximation has a striking consequence 
in how the unitarity (the black disc) limit is approached. In the GM case unitarisation happens because 
of cancellations between contributions of non-saturating multiple gluon exchanges, while in the BK case 
multiple scattering effects are contained in the single gluon density that saturates at a certain small value 
of X. These differences in the mechanism of unitarisation do not affect, however, the crucial qualitative 
feature of the dipole cross-section: the transition from a power-like growth with decreasing x in the 
colour transparency regime to a flat (possibly ~ ln(l/x)) behaviour in the black disc limit. Thus, the 
necessary modelling of the dipole cross section for large dipole sizes is strongly constrained. 

A third type of parameterisation for the dipole cross section does not assume any mechanism for 
unitarisation. It is a two-component Regge model (FS04 Sat) [7], which uses hard Pomeron behaviour 
(~ fQj. sjjiall dipole sizes r < tq and soft Pomeron behaviour (~ x^^'^) for large dipole 

sizes r > ri, with linear interpolation between the two regions. Again, a factorising impact parameter 
dependence is assumed. Saturation effects are modelled by allowing tq to move to lower values with 
decreasing x. This feature was found to be preferred by the HERA data, whereas the two-component 
Regge model with a fixed tq was disfavoured [7]. 

We compare the dipole model parameterisations in Fig. [3^, where the 6-integrated dipole cross 
sections have been divided by the trivial factor in order to emphasise the differences at small r. We 
restrict attention to dipole model parameterisations which have been shown to give a good fit (with charm 
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quarks included) to recent HERA inclusive structure function data, meaning a per data point of 1. 
This excludes, for example, the original GBW parameterisation [4] and the unsaturated two-component 
Regge model [7]. All parameterisations shown in Fig.[3K are similar at intermediate dipole sizes where 
they are most constrained by HERA data. At very small dipole sizes the b-Sat model deviates from 
the other parameterisations, as it is the only one which incorporates explicit DGLAP evolution. The 
b-Sat model was found to be preferred over the b-CGC model for observables sensitive to relatively 
small dipole sizes [12]. There are also differences between the parameterisations in the approach to the 
unitarity limit at large dipole sizes. For example, the b-Sat and b-CGC dipole cross sections tend to a 
constant at large r only for a fixed b, but not when integrating over all impact parameters. 

In order to compare the magnitude of unitarity corrections between various models it is customary 
to define a model-independent saturation scale Q|, that is, the momentum scale at which the dipole- 
proton scattering amplitude M becomes sizable. There is no unique definition of Q| and various choices 
are used in the literature. We define the saturation scale Q'g = 2/r|, where the saturation radius rg is 
the dipole size where the scattering amplitude 

AA(x, r5[, 6]) = 1 - e"^ ~ 0.4, (5) 

chosen to match the corresponding quantity, Qg, in the GBW model [4]. Note that this "saturation scale" 
is still far from the unitarity limit where M = 1. The model-independent saturation scale Q| is shown 
in Fig.[3j5: it is generally less than 0.5 GeV^ in the HERA kinematic regime for the most relevant impact 
parameters b ~ 2-3 GeV^^ [1 1, 12]. It should be remembered, however, that any observable will depend 
on integration over a range of dipole sizes, therefore even at high there will be some contribution 
from large dipole sizes r > rs- Moreover, dipole models incorporating saturation fitted to HERA data 
may be extrapolated to very low x and to predict cross sections for nuclear collisions where the saturation 
scale is enhanced by A^^^ [13]. In these situations, multi-Pomeron exchange may become important and 
extrapolation based on single-Pomeron exchange would be unreliable. 



4 Theory outlook: saturation beyond the BK equation in a statistical picture 

The BK equation describes unitarity corrections in the asymmetric configuration, when the target is 
extended and dense and the projectile is small and dilute. In a more symmetric situation, like 7*((5^)p 
scattering at low Q^, the BK approximation is no longer sufficient. In the diagrammatic formulation, 
besides the fan diagram one should then take into account diagrams with closed Pomeron loops. To 
construct a fully reliable and practical theoretical treatment of this complex case has turned out to be a 
prohibitively difficult task so far. Fortunately, the key properties of solutions of the BK equation in the 
low momentum region follow from its universal features and do not rely on the details of the equation. 

In the Kovchegov derivation of the BK equation [25] one uses the Mueller dipole cascade pic- 
ture [2] of the small x QCD evolution. The equation expressed in terms of the dipole scattering ampli- 
tude, Nuv{y) = AA(x, r, h), with Y = ln(l/x), reads 

^ = ^ [d^w ^""'J^' [N^^ + N^^ - N^^ - N^^N^.] (6) 

oY 271 J [u — w)^{w — v)^ 

where u = h — r/2, and v = b+r /2 (assuming z=l/2 in the definition of b, cf. Sec. 2). The equation has 
two fixed points: the repulsive one, N^v = 0, from which the solution is driven out by the linear term, and 
the attractive one, N^v = 1> where the linear and nonlinear term compensate each other. This scenario 
of linear growth of the amplitude tamed by non-linear rescattering effects is common to all existing 
approaches to the saturation phenomenon. In the uniform case, when N does not depend on the impact 
parameter, b, this combination of growth and nonlinearity was shown to lead to a geometric scaling 
property [28] of the solutions, Nuv{Y) = N{\u — v\'^Ql(Y)) for y S> 1, irrespective of the initial 
conditions [29]. For the 7*p cross section, geometric scaling implies that 
which was observed in HERA data [28]. 
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Interestingly enough, the geometric scaling property of the BK equation does not depend on the 
details of either the linear or the non-linear term. Therefore the scaling is a robust and universal phe- 
nomenon. In particular, the BK equation belongs to the same universality class as a simpler and well un- 
derstood Fisher-Kolmogorov-Petrovsky-Piscounov (FKPP) equation [29], dtu{x, t) = d^^u + u — u^, 
where the rapidity is mapped onto the time t and the logarithm of the dipole size onto the real variable x. 
Employing this connection it was proved that, indeed, both the emergence of geometric scaling and the 
rapidity evolution of the saturation scale are universal phenomena and do not depend on the details of 
the BK equation, provided that the initial condition is uniform in the impact parameter space. 

The statistical framework implied by the Mueller dipole model may also be used to provide some 
quaUtative insight into the effect of "Pomeron loops" in the scattering amplitudes [30, 31]. This effect 
corresponds to a stochastic term added to the FKPP equation [31], 

dtu{x, t) = d^^u + u — + ^/u{l — u)r] (7) 

where r] is the white noise. The origin of stochasticity can be traced back to the discreteness of the colour 
dipoles in the Mueller cascade model. The BK equation is derived in the mean field approximation when 
the density of colour dipoles in the projectile is large enough (n » 1) that statistical fluctuations in 
the number of dipoles can be neglected. In this case, Nuv is an averaged dipole scattering amphtude. 
At the edge of the dense regime of the dipole distribution, however, the dipole occupation number is 
small, n ~ 1, so the statistical fluctuations play an important role. It was reahsed in Ref. [30] and 
subsequently developed in Ref. [31] that these fluctuations get enhanced in the y-evolution and affect 
the global properties of the amplitude. In this approach the saturation scale becomes a stochastic variable 
that fluctuates from one scattering event to another, with a lognormal distribution with the variance 
= DY, where D ~ as/ln^{l/a^) [32]. The most important result of fluctuations is a new scaUng 
of the physical amplitude, called diffusive scaling [31]. Namely, the dipole scattering amplitude Nuv{y), 
should depend only on one variable, ^ = (In(r^) + (in Qiy) / \/ DY . Note that the factor \/DY in the 
denominator which spoils the geometric scaling is of the diffusive origin. A first attempt to trace the 
diffusive scaling in the HERA data on F2 was presented in Ref. [33] with a negative result. This would 
suggest that Pomeron loops introduce only a small effect in the HERA data. 

The results presented here neglect the impact parameter dependence of the scattering amplitudes, 
assuming that the high energy QCD evolution is local in the transverse coordinate space. Thus the local 
evolutions at different 6's are uncorrelated. Recent numerical studies [34] suggest that this is a quite 
accurate picture of high energy scattering if the dipole size is significantly smaller than the target size. 

Recently, an interesting attempt was made [35] to explicitly model the colour dipole cascade taking 
into account effects related to Pomeron loops. In more detail, subleading effects in the 1/Nc expansion 
were phenomenologically incorporated that lead to a possibility of colour dipole reconnections in the 
dipole wave function. The resulting dipole-dipole scattering ampUtudes were shown to respect with 
good accuracy the symmetry between the target and the projectile, which does not hold in the absence 
of the colour reconnection. The approach employs Monte-Carlo methods and was shown to be quite 
successful in describing total cross-sections and many diffractive observables. 

5 Concluding remarks 

The dipole models applied to HERA data on inclusive and diffractive processes provide a successful 
unified description of most observables. These analyses provide significant evidence for sizable unitarity 
(rescattering) corrections to the single scattering approximation, that is used in the linear QCD evolu- 
tion equations, in both DGLAP and BFKL. These corrections become strong below the saturation scale, 
Qs{x). The determination of the saturation scale within different dipole models yields consistently that 
Qs < 1 GeV, over the HERA kinematic range. Qs is found to increase with 1/x, approximately as 
Q'gix) (l/x)'^*' with Xs ^ 0.12-0.2, depending on the model. Both these properties of Qs suggest 
that the onset of perturbative saturation is probed at HERA, and that non-perturbative effects may still be 
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significant around Qs- Fortunately, the key results on the saturation phenomenon obtained within per- 
turbative QCD are universal and should remain vaUd despite a possible non-perturbative contamination. 
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